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ABSTRACT 

We present deep near-infrared (NIR) medium-bandwidth photometry over the wavelength range 
1-1.8/im in the All- wavelength Extended Groth strip International Survey (AEGIS) and Cosmic 
Evolution Survey (COSMOS) fields. The observations were carried out as part of the NEWFIRM 
Medium-Band Survey (NMBS), an NOAO survey program on the Mayall 4m telescope on Kitt Peak 
using the NOAO Extremely Wide-Field Infrared Imager (NEWFIRM). In this paper, we describe the 
full details of the observations, data reduction and photometry for the survey. We also present a public 
if-selectcd photometric catalog, along with accurate photometric redshifts. The rcdshifts arc com- 
puted with 37 (20) filters in the COSMOS (AEGIS) fields, combining the NIR medium-bandwidth 
data with existing ultraviolet (UV; Galaxy Evolution Explorer)^ visible and NIR (Canada- France- 
Hawaii Telescope and Subaru) and mid-IR (Spitzer/IRAC) imaging. We find excellent agreement 
with publicly available spectroscopic redshifts, with a-z/{l + z) 1-2% for ^4000 galaxies at z = 0-3. 
The NMBS catalogs contain ~ 13,000 galaxies at z > 1.5 with accurate photometric redshifts and 
rest-frame colors. Due to the increased spectral resolution obtained with the five NIR medium-band 
filters, the median 68% confidence intervals of the photometric redshifts of both quiescent and star- 
forming galaxies are a factor of ^ 2 times smaller when comparing catalogs with medium-band NIR 
photometry to NIR broadband photometry. We show evidence for a clear bimodal color distribution 
between quiescent and star-forming galaxies that persists to z '--^ 3, a higher rcdshift than has been 
probed so far. 

Subject headings: catalogs - galaxies: distances and redshifts - galaxies: high-redshift - surveys 



; 1. INTRODUCTION 

The growing number of deep, wide-field multiwave- 
I length surveys over multiple bands in the past decade has 
led to new insights into the physical processes that gov- 
ern galaxy formation and evolution. With the advent of 
larger ground-based telescopes with improved NIR cam- 
eras and a new generation of space-based telescopes, sur- 
veys have pushed to higher redshifts and further down 
the luminosity and mass functions. Although broadband 
photometric surveys of statistically significant samples 
of galaxies now span large portions of cosmic time, there 
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exists a fundamental limitation: the accuracy of the pho- 
tometric rcdshift estimates. Accurate photometric red- 
shifts are necessary to measure reliable rest-frame colors, 
stellar population properties, and the environmental den- 
sities. 

Deep optical medium-band photometric surveys 
such as the Classifying Obj ect by Medium -Band 
Observations- 17 (COMBO- 17; IWolf et all [2003) and 
the C osmic Evolution Survey (COSMOS) with 30- 
bands (jllbert et al.ll2009() have photometric redshifts ac- 
curacies of crAz/{i+z,p^a) ^'^'^ large samples of galax- 
ies out to z ~ 1. By using medium-bandwidth filters, 
these surveys have improved the accuracy of photomet- 
ric redshifts by factors of 2-3. 

As prominent rest-frame optical spectral features - 
notably the 4000A break- are shifted into the near- 
infrared at z > 1.5, most studies at high rcdshift have 
focused on those galaxies that are relatively bright at ob- 
served optical (res t -frame ultraviolet) wavelengths (e.g., 
ISteidel et al.l 119961 I1999() . Consequently, t hese studies 
have missed relatively red galax ies (see, e.g.. lFranx et al.l 
[2ooa I van Dokkum et all l2006f ). Concessions must be 
made for studies of complete samples of high redshift 
galaxies, where one ca n either work with small, bright 
spectr oscopic samples (jCimatti et al.l 120021 : iKriek et al.l 
I2008H ). or use large samples that rely on less accurate 
broadband photometric r edshifts (e.g. [Dickinson et al.l 
[2OOI iFontana et aI1[200l and many other studies). 

Here, following previous medium-band surveys in the 
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Fig. 1. — A small region of the COSMOS field (roughly 2' on a side), shown for a broadband J-filter (left) and a three color image from 
the medium-band J1J2J3 filters (right). The red object in the bottom left of the image highlights the spectral features that we are able to 
resolve within only 0.2/im in wavelength (see inset SED). The J-filtcrs trace out the Balmer/4000A break with higher resolution than the 
standard broadband filter, allowing for an accurate photometric redshift of this quiescent galaxy at 2:=2.02. The normalized transmission 
curves of the filters in the NMBS catalogs are shown in the top panel; the additional COSMOS filters (e.g. the Subaru optical medium-band 
filters) are shown in grey-scale behind the five optical filters. 



optical, we extend the redshift range of accurate photo- 
metric rcdshifts out to z ~ 3 using medium-band NIR 
filters. The NEWFIRM Medium-Band Survey (NMBS) 
is a 75 night NOAO survey program on the Kitt Peak 
4m telescope. 

An outline of the paper follows. We introduce the de- 
tails of the NMBS in § [5] and describe the image pro- 
cessing and optimization in § |31 The source detection 
and photometry are next described in § 21 elaborating on 
all of the details involved in extracting the NMBS cata- 
logs. In § I4.10[ we compare the NMBS catalogs to other 
publicly available catalogs within the fields, finding good 
agreement. We derive photometric redshifts and rest- 
frame colors in § [5] and show that the accurate redshifts 
and colors enable us to identify quiescent galaxies out to 
z ^ 3 in § |6l Finally, we demonstrate the improvements 
enabled by the medium-bandwidth filters relative to the 
standard broadband filters by comparing the confidence 
intervals of the photometric redshifts in § [71 A summary 
of the survey can be found in § |S1 

We assume a ACDM cosmology with 177\/=:0.3, 
r2A=0.7, and Hq—IO km Mpc"-'^ throughout the pa- 
per. All magnitudes are given in the AB system. 

2. THE NEWFIRM MEDIUM-BAND SURVEY 

The NMBS employs a new technique of using medium- 
bandwidth NIR filters to sample the Balmer/4000A 
break from 1.5 < z < 3.5 at a hig her resolution than 
the s tandard broadband NIR filters (jvan Dokkum et ahl 
l2009l) . thereby improving the accuracy of photometric 



redshifts. A custom set of five medium bandwidth filters 
in the wavelength range of 1-1.8/im were fabricated for 
the NEWFIRM camera on the Mayall 4m telescope on 
Kitt Peak for the NMBS. The Ji-band is similar to the 
y-band, the canonical J-band is split into two filters J2 
and J3, and the i?-band is split into two filters Hi and 
H2 (see the top panel in Figure H]). T he full details of the 
mediu m-band filters can be found in Ivan Dokkum et al.l 
(12001). 

Figure [T] demonstrates the power of the medium-band 
filters, comparing a traditional J-band image to a three- 
color image comprised of the J1J2J3 filters. Although 
the wavelength range of the three medium bandwidth J- 
filters only covers 0.2/im, the filters are able to resolve 
strong spectral features such as Balmer/4000A breaks or 
emission lines. For example, J1J2J3 trace out the 4GG0A 
break of the massive, quiescent galaxies at z = 2.02 
shown in Figure [1] When using the broadband J and 
H filters alone, the photometric redshift uncertainty is 
in (zNMBS - ^broadband)/ (1 + ZNMBs), typical of the 
inherent uncertainties associated with broadband photo- 
metric redshifts. Through sampling the Balmer/4000A 
break region of the spectral energy distributions (SED) of 
the galaxies at 1.5 < z < 3.5 with higher resolution, the 
uncertainties in the photometric redshifts and rest-frame 
colors decreases by about a factor of two (sec §|7]). 

2.1. Field Selection 

The survey targets two fiel ds wit hin the COS- 
MOS (;Scoville..2007;) and AEGIS (jDavis et al.,,2007.) sur- 
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veys, chosen to take advantage of the wealth of publicly- 
available ancillary data over a broad wavelength range 
(see §[3H). The NEWFIRM 27'.6x27'.6 pointing within 
the COSMOS field is centered at a = 9''59™53.3", S = 
+02°24'"08^ (J2000). This point ing overlaps with the 
zCOSMOS deep redshift survey fjUWy et al.l[2007h and 
the upcoming Ultra- VISTA survejO. The pointing in 
the AEGIS field is centered at a 14''18'"00", S = 
-K52°36™07^ (J2000). 

2.2. Observations 

The observations were carried out on the Mayall 4m 
telescope on Kitt Peak using the NEWFIRM camera 
with the J1J2J3H1H2 medium bandwidth filters and the 
K broadband filter. Data were taken over the 2008A, 
2008B and 2009A semesters, for a total of 75 nights. The 
two main fields were observed whenever conditions were 
reasonable (i.e., no significant cirrus and seeing typically 
<1.5") and they were available at an airmass <2 and 
>20° away from the moon. 

Due to the dominance and short-term variability of 
the sky background, ground-based NIR imaging requires 
many short dithered exposures. T his observing pro- 
cedure is w ell established (see, e.g., iLabbe et al.l 120031 : 
iQuadri eFa l. 2007), and facilitates the removal of arti- 
facts remaining from bright objects that can leave resid- 
ual images in subsequent exposures. The telescope fol- 
lows a semi-random dither pattern with a dither box of 
90", enabling good background subtraction without sig- 
nificantly reducing the area with maximal exposure. The 
dither box is sufficiently large to "fill in" the gaps be- 
tween the four NEWFIRM arrays. 

At each dither position, the exposure times were typi- 
cally 1 X 60 seconds (coadds x individual exposure time) 
for J1J2J3, 3x20 seconds for H1H2 and 5x12 seconds for 
K. The internally-coadded frames make longer exposure 
times possible at each dither position while staying in 
the linear regime of the array. 

We obtained a total on-sky integration time between 
17 and 31 hours in each filter in each field over the 3 
semesters (see Table [1] for details). This represents all 
usable data, including aborted sequences or frames that 
were outside weather specifications or seeing constraints; 
the total time on-sky was 294 hours. When tallying only 
those frames that get included in the final combined im- 
ages, the total integration times range from 12 to 24 
hours (a loss of ^10-30% of the data). 

An overview of the median seeing and background lev- 
els measured in each observing sequence (ranging from 
10-90 minutes, typically 60 minutes) is given in Figure[2l 
The median seeing is measured from bright stars in the 
individual frames that are fit with a two-dimensional 
Moffat model using the robust, non-linear least squares 
curve fitting algorthm mpf it in IDL (Markwardt 2009). 
The data are of excellent quality with clean images that 
have no obvious systematic problems. 



3. DATA REDUCTION 
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Fig. 2. — The median seeing, as implied by a Moffat fit, and sky 
baekground in the combined sequences (typically ~60 minutes) for 
Ji , J2 and J3 (light grey) , Hi and H2 (grey) , and K (black) . 

The data were reduced using a custom IDL code writ- 
ten by one of us (IL), following the method of the IRAjQ 
package XDIMSUMg] . The re duction process largely 
follow s ILabbe et al.l (I2003D . iForster Schreiber et al.l 
(l2006l) . IQuadri et all (I2007D and lTavlor et al.l (|2009bD . 
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Below, we summarize the process, highlighting modifi- 
cations to the standard process described in the papers 
listed above. 

3.1. Image Processing 

Every individual raw frame was visually inspected to 
identify any poor quality frames. To fall in this catagory, 
the frame must either (1) contain artifacts such as satel- 
lite streaks, ghost pupils or other transient features on 
the detector, (2) have a strongly varying background 
across the detector, (3) have severe elongation of the 
point spread function (PSF) or (4) have double images 
of all objects due to "jumps" of the telescope or wind 
shake. Between 10-30% of the raw data was discarded 
for each filter, resulting in final combined mosaics with 
negligible artifacts. 

The dark current is nonneglibible in NIR cameras and 
thus dark images with the appropriate exposure times 
relative to the co-added science images are subtracted. 
Next, the raw images must be divided by the fiat field 
image to remove variations in the sensitivity across the 
array. 

Wc constructed fiat fields from dome flats in combina- 
tion with observations of the open cluster M67, account- 
ing for both the small and large scale variations across the 
arrays, respectively. First, we determined the pixel-to- 
pixel response from dome flats. The large scale structure 
in the dome flats was removed by median flltering using 

^■^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

Ihttp: //iraf . noao . edu/iraf /ftp/ftp/extern/ xdlmsum/ 
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TABLE 1 

Summary of Observations 



Field 


Filter On-Sky Tirnc''[lir] 


Integration Time'^[hr] 


Image Quality^ [ 


] 5(T Depth^ Aperture Correction*^ Zero Point 




Jl 


29.3 


25.0 




1 1 Q 


25.2 


1.9 


23.31 




J2 


25.9 


24.1 




l.io 


25.3 


1.9 


23.35 




J3 


29.1 


22.9 




1.08 


24.5 


1.8 


23.37 




HI 


26.5 


22.3 




1.10 


24.1 


2.0 


23.59 




H2 


20.8 


17.1 




1.06 


24.4 


1.8 


23.61 




K 


22.2 


12.8 




1.08 


24.2 


1.9 


23.85 


COSMOS 


Jl 


31.2 


24.9 




1.19 


25.1 


2.0 


23.31 




J2 


21.9 


19.6 




1.17 


24.8 


1.9 


23.35 




J3 


27.5 


24.8 




1.12 


24.7 


1.9 


23.37 




HI 


22.3 


14.7 




1.03 


24.2 


1.8 


23.59 




H2 


20.0 


17.0 




1.24 


24.0 


2.1 


23.61 




K 


17.0 


11.8 




1.08 


24.2 


1.9 


23.85 



" The total on-sky time represents all usable data, including aborted sequences and frames outside of the weather specifications or seeing 
constraints. 

^ The total integration time represents only those frames that were included in the final combined mosaics. 

The direct measurement of the FWHM from IRAF, independent of a profile model. The derivative of the enclosed flux profile is computed 

and the peak of the azimuthally averaged radial profile is found, where the FWHM is twice the radius of the profile at half the peak value. 

These are 5a depths in a circular aperture of 1.5" diameter, which corresponds to a total magnitude that is about 0.7 mag brighter for point 

sources for a seeing of 1.1". 
" The aperture correction for a point-source to convert the 5cr depths to total magnitudes. 



a 51x51 pixel box, fitting a 7th order 2D Legendre poly- 
nomial to the median filtered image, and dividing by this 
fit. It was verified that subtracting rather than dividing 
by the fit changes the resulting residual image by <1% 
(except for the iV'-band, where differences can be up to 
^--^5% in the corners of the arrays). 

The simple flat fields are not enough to correct to the 
large-scale sensitivity of the array due to imperfect illu- 
mination corrections. We therefore determine the sensi- 
tivity variations across each array from observations of 
the open cluster M67 in an 8 x 8 grid of positions spanning 
28'x28'; this way, many stars were observed on many ar- 
ray positions. The M67 data were dark subtracted, flat 
fielded using the "fiattened" dome fiat (i.e., the dome fiat 
with large scale structure removed) and sky subtracted. 
The fiuxes of stars in M67 were m easured using Source 
Extractor (jBertin fc Arnoutill996D . with a 4" diameter 
aperture. Using the RA and DEC of the detected objects, 
sequences of unique stars imaged on multiple positions 
were identified. From these stars, a response image was 
generated refiecting the large scale response of the arrays. 
The scatter in the best-fit response was typically 1-3%, 
with larger errors possible in small areas for some ar- 
ray/filter combintations. These "response images" were 
then multiplied by the fiattened domefiat to create the 
fiat fields. 

The individual dark-subtracted and fiat-fielded frames 
are combined in two passes. During the first pass, the 
sky background is created for each science frame from 
a running median of the dithered sequence of the four 
preceeding and four subsequent science images. Next, 
the dithered sky-subtracted images are cross-correlated 
with the central frame to determine the relative shifts of 
the individual frames. The relative flux scaling between 
the images is determined by measuring the flux of stars 
spread across the field. Finally, the images are shifted 
to a common reference frame using linear interpolation 
and are combined. All objects in the combined image are 
detected using a simple thresholding algorithm to gener- 
ate an object mask. This mask is then shifted back into 
the reference frame of each individual science exposure 



and the background subtraction process described above 
is repeated in a second pass, using a mean combine with 
cubic interpolation. The second pass is optimal for the 
background subtraction as the objects are masked out 
during the calculation of the running median. 

3.2. Further Optimization 

In the reduction, we mask out any dead or hot pix- 
els (or generically "bad" pixels) using a custom master 
bad pixel file. Bad pixels have a non-linear response and 
these bad pixels can be easily separated from "normal" 
functioning pixels during twilight. The signal of good 
pixels increases or decreases linearly with time, and so 
we fit the average fiux of each quadrant as a function of 
the time with a linear function. Any dead or hot pix- 
els will exhibit variable behavior, resulting in drastically 
different best-fit slopes compared to properly function 
pixels. We fiag any pixels with deviations >15%. The 
final bad pixel file contains a flag for any pixel that has 
been flagged in at least two of the twilight sequences, 
in addition to the original bad pixel file from the NEW- 
FIRM calibration library. Although we are conservative 
with the bad pixel file, the total fraction of bad pixels is 
still small (^ 7%). Any remaining bad pixels in each in- 
dividual science frame are removed with a sigma-clipping 
algorithm before combining the images within the reduc- 
tion. 

The nonlinearity of the NEWFIRM arrays increases 
with additional signal, reaching ^--^5% at 10,000 ADU. 
Beyond this level, the nonlinearity increases very steeply 
as the detector approaches saturation at 12,000 ADU. 
Coadding frames allows us to keep exposure times short 
thereby avoiding this nonlinear regime, as does observing 
when the sky levels are relatively low in a particular fil- 
ter. A nonlinearity correction is applied to all individual 
frames, using the following empirical equation: 



where x is the non-linear (raw) frame, y is the non- 
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hnearity corrected frame, n is the number of internally 
coadded frames and 6.64 x 10~^ is the average coefficient 
for the four arrayQ. The non-linearity correction to the 
photometry is typically < 1% for the majority of frames. 

Another issue in NIR detectors is persistence images 
of bright objects, where the array has not completely 
recovered from the previous exposures. To alleviate this 
problem, we create a second object mask of the cores 
of the brightest objects. These are used to mask the 
pixels that contained the brightest objects in the previous 
two exposures in the final pass of the reduction. The 
detectors also exhibit bias residuals which are constant 
along rows. After masking all sources, the residual bias 
is removed by subtracting the median of each row in the 
individual sky-subtracted frames. 

The final major optimization of the image combination 
includes weighting the individual sky-subtracted frames 
to maximi ze the signal-to-noisc (S/N) ratio (sec e.g., 
iLabbe et al . 2003; Quadri ct al. 2007). This method sub- 
stantially improves the final image depth and quality by 
assigning weights to the individual frames that take into 
account variations in the seeing, sky transparency, back- 
ground noise and PSF ellipticity. To optimize the S/N, 
the weight is proportional to the square of the flux scale 



|http: //wMw .noao ■ edu/staf f /med/newf Irm/f ewompt/ 
|newf irm_lln6arlty . pdf] 



(scale-;) divided by the median sky background (sky^) and 
the median size of the seeing disk (FWEM^), penalized 
for ellipticity: 

scale^ 

~; FWHM^ 

3.3. Astrometry 

The weighted, sky-subtracted individual frames are 
combined into a single mosaic with a pixel scale of 0.3" 
pix^^ (resampled from the native NEWFIRM pixel scale 
of 0.4" pix~^). The combined NMBS mosaics arc aligned 
with the CFHT /-band images (see ij3.6p . with an astro- 
metric precision of ^0.1" over the entire field of view. 
The relative registration between the NMBS and CFHT 
images must be precise to allow for accurate colors and 
cross-identification of sources. 

The initial registration within the code is performed 
using a linear transformation in the first pass of the re- 
duction. The rms variation in the position of individual 
sources is about 0.1-0.2" (0.3-0.7 pixels) in most cases, 
but can be as high as >0.5-l" (2-3 pixels) in the corners 
of the arrays. We fit the residual distortions with a sec- 
ond order polynomial, and combine this fit with the origi- 
nal transformation. This combined distortion correction 
is then applied to the images, which are registered us- 
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Fig. 4. — The AEGIS X-band image (roughly ~ 30 X 30 arcminutes), with the inset demonstrating the quality and depth of the data. 

3.5. Noise Properties 

The final combined Ji, J2, J3, Hi, H2 and K images 
constructed from tfie individually registered, distortion- 
corrected, weighted averages of all frames are of excellent 
quality. The flatness of the background is readily visible 
in the final combined /f-band images shown in Figures [3] 
and|4l Furthermore, the relatively deep and wide images 
arc rich in compact sources. The combined images are 
slightly shallower near the edges and center of the image 
where there are gaps between the NEWFIRM arrays, 
as these areas receive less exposure time in the dither 
pattern. This is reflected in the weight maps containing 
the total exposure time per pixel. The reduced NMBS 
images and weight maps will be made publicly available 
through the NOAO archivcEl. 

Well-understood noise properties are important both 
for calculating the depths of the final combined images, 
as well as measuring accurate photometric uncertainties. 
Although the noise properties of the raw images are well 
described by the variance of the signal per pixel, the 
reduction process introduces correlations between neigh- 
boring pixels. Additionally, small errors in the back- 
ground subtraction may contribute to the noise. 

To estimate the noise, we follow iLabbe et all (|2005f ) 
and empirically fit the dependence of the normalized me- 



|http : //archive ■ noao . edu/nsaZ 



ing a cubic interpolation. The resulting rms variation in 
position of the individual sources in the final combined 
images is <0.1" (0.3 pixels). 

3.4. Photometric Calibration 

NIR spectro-photomctric standard stars from the Cal- 
ibration Database System (CALSPECjB were observed 
in wide five-point dither patterns on photometric nights 
in the 2008A semester. All of the other data were directly 
tied to the photometric data from 2008A. Synthetic mag- 
nitudes of these stars were calculated by integrating their 
observed (HST/NICMOS) spectra ii i the five medium- 
band width filters (see Table 2 in Ivan Dokkum et aD 
l2009f ). We adopt the mean zero points from these ob- 
servations, listed in Table [TJ The zero points derived 
from these standard stars arc remarkably stable, with 
variations of <0.02 magnitudes throughout the observ- 
ing program. 

A Galactic extinction correction has been applied to 
the photometry with in the catalogs, as estimated from 
iSchlegel et all (|1 9981) . The corrections ranged from 4.5% 
to 0.3% for the u-K photometry in the AEGIS field, with 
slightly larger values of 8% to 0.6% in the COSMOS field 
due to its lower Galactic latitude. 



^ |http: //www. stsci . edu/hst/observatory/cdbs/calspec ■ 
I html 
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dian absolute deviation (nmad) of the background as a 
function of the linear size N of an aperture with area A, 
where N = \fA. We measure the fluxes in >1000 inde- 
pendent circular apertures randomly placed on the reg- 
istered, sky-subtracted, convolved images (see § 14.11 for 
a description of the PSF matching) in a range of aper- 
ture sizes. Any apertures that overlap with the detection 
image segmentati on map are excluded . 

As described in lQuadri et al.l (j2007[) . the scaling of the 
background noise with aperture size will be proportional 
to iV^ in the limiting case of perfect correlations between 
the pixels within an aperture. In the other extreme case 
of uncorrelated adjacent pixels within an aperture, the 
background noise scales with the linear size N . The latter 
standard relation is unrealistic given the correlations be- 
tween pixels that were introduced by the imperfect back- 
ground subtraction, but also by undetected sources, re- 
sampling of the pixels, etc. The true scaling between the 
background and the linear size falls somewhere bet ween 
these two extremes. Following iQuadri et al.l ()2007D . we 
parameterize the noise properties as: 

O-nmad = diaiV^, (3) 

where 1 < /3 < 2, (Ti is the standard deviation of the 
background pixels and a is the normalization. In Fig- 
ure El we see that the width of the noise distribution 
within an aperture increases with increasing aperture 
size. The best fit scaling of the background fluctuations 
falls roughly between the two extremes of no correlation 
and a perfect correlation from pixel to pixel, with a best 
fit power-law index of /3 = 1.37 in COSMOS (/3 = 1.35 
in AEGIS). 

The background noise measurements are used to esti- 
mate the depths of the images. In Table [U we list the 5(t 
limiting depths. These values are calculated using the 
background fluctuations within a 1.5" "color" aperture 
(see §|4T]). An additional aperture correction determined 
from the growth curves must be applied to these values 
to convert to total fluxes (see Column 7 of Table [1]) . The 
typical aperture correction for a point-source is 1.9 for a 
typical seeing of ~ 1.1". 

3.6. Ancillary Images 

The survey catalog combines the data taken with 
the NEWFIRM camera with publicly-available optical 
ugriz broadband images of both survey fields from the 
Deep Canada-France-Haweaii Telescope Legacy Survey 
(CFHTLS£3L 

using imag e versions reduced by the CARS 
team ()Erben et all 120091 : iHildebrandt et al.l 120091 ) . and 
JHKs broadband imaging from the WIRCam Deep Sur- 
vey (WIRDS; Bielby et al.,m preparation). Additionally, 
we include deep Subaru images with the BjVjr^i'^ z'^ 
broadband filter s and 12 medium-ba nd filters in the 
COSMOS field (iTaniguchi et al.l [20071 Y. Taniguchi et 
al. 2008, in vrevarationS^. We include mid-IR im- 
ages in the Spitzer-lYlAG bands over the entire COS- 
MOS pointing_that are pr ovided by the S-COSMOS sur- 
vey () Sanders et al.l [2007f l. and partial coverage of the 
AEGIS pointing 0. 15 deg^) overlapping with the Ex- 
tended Groth Strip (jBarmbv et al.l l2008| i. The cata- 

1'^ |http : //mmm. cf ht . hawali . e du/Science/CFH TLS7] 

http : //irsa . ipac . caltech . edu/data/COS MOS/lmagesT] 



0.7 1 ' ' ' I ' ' ' I ' ' ' I " ' ' I ' ' ' I ' ' ' I ' ' ' I ' 




0.0 r ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ ■ ■ I ■ 
2 4 6 8 10 12 14 

N [pixels] 



Fig. 5. — The nmad value of background fluctuations for the 
COSMOS K-hand image within an aperture with area A as a func- 
tion of the Hncar size N , where N = \fA. The filled triangles il- 
lustrate the measured values, whereas the black solid curve is the 
best fit to the data following the parametertization in Equation |3] 
The grey solid curve is the best flt to the AEGIS data (not shown 
here). The dashed curves show the expected scaling relations of 
the background fluctuations in the case of no correlations amongst 
the pixels (bottom) and perfect correlations of all pixels within the 
aperture (top). Similar fits were performed for all fllters in both 
fields. 

logs include 5piteer-MIPS fluxes at 24/im that are de- 
rived from images provided by the S-COSMOS and FI- 
DElF^ surveys, with the same coverage as the IRAC 
data. Finally, we incorpora te NUV and FUV Galaxy 
Evolution Explorer (GALEX: rMartin et aT]|2005D data in 
both flelds from the Multimission Archive at Space Tele- 
scop e Science Institute ( MAST) into the catalogs (see 
also iZamojski et"ani2007t for a description of the COS- 
MOS GALEX data). 

All ancillary images arc matched to the same pixel 
scale (0.3" pix'^) and pointing of the NMBS images us- 
ing the IRAF task wregister. The details of the source 
detection and extraction of photometry can be found in 
§ 14.21 In general, for objects that have flat spectra in 
F\, all of the ancillary data are deeper than the NMBS 
images with typical FWHM values of 0.8" for the opti- 
cal images. For consistent photometry in all (optical to 
NIR) bands, the images have all been convolved to the 
broadest PSF as determined from the growth curves (see 
^ 14.11 for more details). 

4. PHOTOMETRY AND SOURCE DETECTION 

We have 20 independent photometric filters in AEGIS, 
and in COSMOS we have 37 (FUV-8 urn). The optical 
and NIR images were convolved to the same PSF be- 
fore performing aperture photometry, so as to limit any 
bandpass-dependent effects. The source detection was 
don e on the PSF matched images with Source Extrac- 
tor (|Bertin fc Arnout^ll996l ) in relatively small apertures 



|http : //irsa ■ ipac ■ caltech . edu/data/SPITZER/FIDEL/ 
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Fig. 6. — The ratio of the growth curve for all optical to NIR bands Fi relative to the band with the worst seeing Fx (where x = Hi in 
AEGIS and H2 in COSMOS). The top panel shows the ratio before convolution (dashed lines) and the bottom panel shows after convolution 
(solid lines). Note the different y-axis scales of the top and bottom panels. The PSF matched images have accurate colors to < 2% for the 
catalog aperture diameter of 1.5". 

ratio should be unity at all aperture diameters. At the 
aperture size used for the flux measurements (solid black 
line), the PSFs of all bands are matched within <2% for 
both fields. 

Given the noise properties of the images (see § 13. 5p and 
the growth curve of the PSFs, we estimate the aperture 
size that optimizes the S/N for point sources in Figure[7l 
This optimization enables accurate color measurements 
which are necessary for photometric redshift calculations 
and modeling of the stellar populations. We adopt an 
aperture of 5 pixels (1.5"), defining the color aperture 
as: 



chosen to optimize the S/N (see ? 13. 5p . To extract the 
photometry from the ultraviolet and mid-IR images, we 
use an alternate source fitting algorithm especially suited 
for heavily confused images for which a higher resolution 
prior (in this case, the /C-band image) is available. 

4.1. Optical and NIR PSF Matching 

The FWHM of the PSF varies significantly between 
the optical and NIR bands: point sources in the opti- 
cal images typically have FWHM of approximately 0.8", 
whereas the point sources in the medium-band NIR im- 
ages have FWHM of 1.0"-1.2". Directly performing 
aperture photometry on these images would lead to dif- 
ferent fractions of the flux falling within the aperture 
for different bands, thereby resulting in systematic color 
errors. 

To minimize this problem, we degrade all of the optical 
and NIR images for each field to the broadest PSF. We 
produce high S/N PSFs for each band from bright (un- 
saturated) stars covering the entire field. Square postage 
stamps (32 pixels on a side, or 9.6") of the PSF stars were 
created, masking any nearby objects and averaged to cre- 
ate empirical PSFs. A kernel is then generated using the 
Lucy-Richardson method algorithm (lucy in IRAF) to 
convolve the empirical PSF to the broadest PSF {Hi in 
AEGIS and H2 in COSMOS). 

Growth curves were determined for the empirical PSFs 
in each band, measuring the enclosed flux as a function 
of the aperture size. In the top panels of Figure [6l we 
show the ratio of the growth curves of each band rela- 
tive to the broadest PSF in both fields before degrading 
the PSF (dashed lines). Generally, the broadest PSFs 
(the slowest-growing growth curves) are the longer wave- 
length bands. The bottom panels of Figure [5] show the 
ratio of the growth curves once they have been smoothed 
to the worst seeing. For "perfect" PSF matching, this 



color, ^ 



(4) 



<0.75" 



where Fr, 



)lor, 



is the fiux for each band i within a circular 
^ of 0.75". The color 



aperture with a radius r = yx^ 

aperture defined in Equation |4] is small enough to give 
nearly optimal S/N for point sources and large enough to 
minimize possible systematics from residual astrometric 
or PSF matching uncertainties. We note that the K- 
band has a relatively flat slope in S/N, whereas the im- 
provements are most prominent for the Ji band (shown 
in Figure [Zl). There is a 20-30% loss of S/N when us- 
ing larger apertures of 2-3" in bands blueward of the 
A'- band. 

The empirical estimate of the background noise cTnmad 
described in § 13.51 was repeated for all optical and NIR 
images within the color aperture of 1.5". From these 
values, we estimate the flux uncertainies for each object 
in each band as: 
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color, ^ 



F, 



color, 2 



^Wi/Wmcd 
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Fig. 7. — The signal-to-noisc as a function of the aperture size, 
derived from the ratio of the growth curve and the noise proper- 
ties for the K (black), J\ (blue) and X-band images (red). The 
solid lines are for the COSMOS field and the dashed lines are for 
AEGIS, with an "x" marking the peak in S/N. The FWHM of the 
PSF matched _fC-band image is shown with an arrow. Wc use an 
aperture diameter of 5 pixels (1.5") when creating the NMBS cat- 
alogs so as to optimize the S/N while remaining larger than the 
FWHM. 

The photometric error for each band i is the sum in 
quadrature of the background uncertainty (cTnmad.i; de- 
fined in Equation [3|) weighted by the fractional expo- 
sure Wj at each object's location relative to the median 
Wmedian,i, and the sourcc Poisson noise Fijgi for each 
band, with the former dominating the noise for faint 
sources and the latter for bright sources. Fi is the flux 
of the object in ADU within a 1.5" aperture and gi is 
the total effective gain. Note that these are the uncer- 
tainties of the aperture fluxes, not the total fluxes. The 
uncertainty in the total fluxes listed in the catalogs is: 



f^total.: 



'^colo 



total 



)lor 



(6) 



The ii'-band total magnitude is calculated from the 
SExtractor AUTO photometry (sec Equation |9]) , with 
an uncertainty calculated as follows: 



'Ktot 



median, 



-?^AUTO 
9K 



(7) 



is the circularized Kron radius (|Kronl 
(Ti, a and j3 are defined in Equa- 



where i?Kion 

floM ) and a I, a and p are dcnnea m 
tion 121 weighted by the fractional exposure 
time. The circularized Kron radius is calcu- 
lated from the semimajor and scmiminor axes as 

i?Kron = \/A.IMAGE X BJMAGE x KRON_RADIUS^ 
using the output SExtractor parameters. The following 
equation can be used to determine the full error in the 
total magnitude for a single band in isolation, including 
systematics: 



^2 

^full.j 



''total. i + '''Ktot "'total, K 



(8) 



The derivation of the uncertainties for the ultraviolet 
and mid-IR data are described in § 14.31 and 



4.2. Source Detection 

The sources are detected in a noise-equalized i^-band 
image using SExtractor in dual-image mode, where the 
aperture photometry is simultaneously measured on the 
PSF matched u-K images. In both fields, the i^-band 
detection image is the unsmoothed (pre-PSF matching) 
NEWFIRM 7^-band image multiplied by the square root 
of the weight map (based on fractional exposure time 
per pixel). The input SExtractor parameters were op- 
timized so as to find all faint objects while limiting the 
number of spurious detections. The optimal setting for 
the minimum number of pixels above the threshold to 
trigger a detection was 5 pixels, with a detection thresh- 
old of 1.2(7 and 32 deblending sub-thresholds. We use 
a Gaussian filter for the detections with a FWHM of 4 
pixels, first performing a global background subraction 
in all optical-NIR images with a background mesh size 
of 32 pixels and a smoothing filter of 3 pixels. We note 
that we used a minimum Kron radius of 3 pixels, com- 
pared to the default value of 3.5 pixels. A fixed aperture 
diameter of 1.5" (5 pixels) was used to both enclose a sig- 
nificant fraction of a given object's fiux while minimizing 
uncertainties due to background fluctuations (see § 14. 1[ 
Figure III). 

Total fluxes were determined using the SExtractor 
AUTO photometry, which uses a flexible elliptical aper- 
ture. This aperture contains the majority of the flux for 
bright objects but may be missing a large fraction of the 
light for fainter objects. Therefore an aperture correc- 
tion is ap plied to convert th e AUTO flux to total flux 
(see, e.g.. iLabbe et al.ll2003D . The Jf-band total flux is 
calculated as follows. 



^''total = ^AUTO X 



R<AUTO 



(9) 



where -R'auto is the AUTO flux and / is the fraction of 
light enclosed within the X-band growth curve generated 
during the PSF matching described in § l4.1l for a circular 
aperture with the same area as the AUTO aperture. We 
force a minimum size of the AUTO aperture equal the 
size of the color aperture, affecting < 0.1% of the sources. 
The aperture correction to the AUTO fluxes is typically 
on the order of 5-10%, extending to higher values for the 
faint, extended sources. The total flux for each band i is 
calculated as follows: 



-^total.2 — -^color.i ^ 



^^total 
color 



(10) 



where Fcoior,i is the flux within a color aperture of 1.5" 
diameter (see Equation 2]) . 

4.3. IRAC/MIPS PSF Matching and Photometry 

The FWHM of the IRAC and MIPS images are a fac- 
tor of ^2 and ^5 broader than the PSF-matched images. 
Furthermore, the IRAC and MIPS images contain sig- 
nificant non-Gaussian structure that results from point 
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source diffraction. Smoothing the optical to NIR pho- 
tometry to the IRAC PSF shape would substantially re- 
duce the S /N and cause significant blending issues in the 
deep optical data in addition to the IRAC/MIPS images. 

We use a source-fitting algorithm developed by one 
of us (IL) that is designed for heavily confused images 
to extract the p hotometry from the IRAC and MIPS 
images (see, e.g. . iLabbe et al.l [20061: l^vts et al.l [20071: 
iMarchesini et al.ll2009l:lWimams et al.ll2009[) . ATsources 
that are bright in K are also typically bright in the IRAC 
images (and to a lesser extent MIPS images), the i^-band 
can be used as a high resolution template to deblend the 
IRAC/MIPS photometry. The information on position 
and extent of the sources based on the higher resolution 
iiT-band segmentation map is therefore used to model the 
lower resolution IRAC 3.6-8.0 /xm and MIPS 24 / xm im - 
agcs. This technique is similar to iLaidler et a l. (2007*), 
with one important difference. ILaidler et al.l ((2007.) use 
the best-fit fluxes directly to calculate flux ratios (col- 
ors), whereas our method adds back in the residuals of 
the fit before performing aperture photometry. 

Each source was extracted separately from the /C-band 
image and, under the assumption of negligible morpho- 
logical K-corrections, convolved to the IRAC/MIPS res- 
olution using local kernel coefficients. The convolution 
kernels were constructed using bright, isolated, unsatu- 
rated sources in the K and IRAC/MIPS bands, derived 
by fitting a series of Gaussian-weighted Hcrmitc func- 
tions to the Fourier transform of the sources. Outly- 
ing or poorly-fit kernels are rejected and a smoothed 2D 
map of the kernel coefficients was stored. All sources 
in each IRAC/MIPS image are fit simultaneously, with 
the flux left as a free parameter. The modeled light 
of neighboring objects is subtracted, thereby leaving a 
clean IRAC/MIPS image to perform aperture photome- 
try in fixed 3" apertures for IRAC and 6" apertures for 
M IPS. For an i ll ustra tion of this technique, see Figure 1 
in lWuvts etldl (|2007D . 

In order to compute a consistent X-IRAC color, we 
measure the flux of the objects on a cleaned iC-band 
image convolved to the IRAC resolution within the same 
aperture to correct the photometry for flux that falls out- 
side of the aperture due to the large PSF size. We then 
scale the photometry to total fluxes by measuring the 
ratio of the total i^T-band flux to the convolved /C-band 
flux within a 3" aperture for each object; 



-FiRAC, catalog = ^IRAC,clcancd(3") X — — — ■ (11) 

-'»-conv('^ } 

where -Firac, cleaned is the flux of the object measured 
with the modeled light of the neighbors removed, A'conv 
is the flux of the object measured in the K-hand image 
after smoothing it to the IRAC resolution, and i^totai is 
the total i^-band flux. The i^iRAC catalog flux is the total 
flux that would have been measured in a 1.5" aperture 
with an aperture correction if IRAC had the same PSF 
as the other data. 

The uncertainties in the measured fluxes are derived 
from the residual contamination of the subtracted neigh- 
bors. The normalized convolved A'-band image of each 
neighboring source is scaled by the formal la error in 
the fitted flux as computed from the covariance matrix 



produced by least-squares minimization. By performing 
aperture photometry on the image of the residual neigh- 
bor contamination, we get the error corresponding to the 
source flux measured within an aperture of the same size. 

The catalog includes the 24 jim flux within flxed 6" 
apertures and total fluxes. To measure the total 24 
/im fluxes, we create a MIPS growth curve from several 
bright, isolated, unsaturated point sources within each 
field. These square postage stamps are 12.6" on a side, 
and we derive an aperture correction from 6" to 12.6" 
of a factor of 2.8. To convert to total flux, we include 
an additional aperture correction for the 17% of the flux 
that falls outside 13'H 

4.4. GALEX PSF Matching and Photometry 

The GALEX Far UV band (1350-1750A) and Near UV 
band (1750"2800A) have 4.5"-6" resolution (FWHM), 
roughly ^--^4-5 times broader than the PSF matched op- 
tical to NIR photometry. The same method described in 
§ 14.31 is used to extract the GALEX photometry. In this 
case, however, the higher resolution prior is the CFHT 
M-band image. 

As the u-band contains many blue galaxies not in the 
A'-band detection image, a new segmentation map was 
generated with SExtractor for the it-band and matched 
to the A'-selected catalog. This allows for all of the ob- 
jects within the field of view to be modeled to clean the 
GALEX images for aperture photometry, regardless of 
whether these objects are in the NMBS A'-selected cat- 
alog. 

Following § 14. 3[ we measure an aperture correction 
from the flux of the objects on the cleaned w-band image 
convolved to the GALEX resolution in order to compute 
a consistent GALEX-m color. We then scale the pho- 
tometry to the same color apertures used for the optical 
to NIR photometry by measuring the ratio of the con- 
volved w-band flux within a 6" aperture and our catalog 
aperture size for each object, analogous to Equation [TT] 

4.5. Deblending 

A signiflcant fraction of the objects detected by SEx- 
tractor are actually blended objects, although classified 
as a single object. Although SExtractor implements a 
deblending algorithm, there is no set of parameters that 
will enable the splitting of close pairs while also detecting 
faint objects, which requires pre-smoothing. Deblending 
the photometry is particularly important when study- 
ing the most luminous galaxies; the fraction of sources 
that are deblcnded increases from ~ 10% at K=22 mag 
to ^ 40% at K=20 mag, reaching as high as ~ 80% of 
sources with K = 18 mag. In this section, we describe 
the method used to deblend the photometry of the single 
"parent" object into its co nstituent "children" , based on 
the technique described bv lLupton et al.l ()2001[ ). 

The deblending is done using the WIRDS Ks-hand im- 
age (0.7" seeing) in COSMOS and a combination of the 
WIRDS ATg-band image (0.7" seeing) where the expo- 
sure is greater than 30% of the maximum and the non 
PSF matched AT-band image in AEGIS (1.1") in the re- 
maining area. Using this criterion, roughly half of the 
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Fig. 8. — An example of the SED of a galaxy that is blended in the original catalog (left) compared to the SEDs of the dcblcnded children 
(right). The inset panel shows the CFHT Ks-hand image with the original and dcblcnded segmentation maps below. The SEDs of the 2 
children are very different. The photometric rcdshift of the first child is usually very similar to that of the parent. 



AEGIS dcblcnding image is from the higher resolution 
WIRDS data. Peaks are found after applying a Differ- 
ence of Gaussians band pass filter that preserves scales 
similar to the seeing size. Object templates are cre- 
ated by comparing pairs of pixels symmetrically around 
the peak and replacing both by the lower value. False 
peaks are rejected if the total template S/N< 3 or if the 
template peak value is less than 2x the combined pixel 
value of all other templates at that location. Finally, 
we solve for the template weights by a linear fit of the 
templates to the image, and then assign flux to every 
pixel of each child according to the relative contribution 
of the weighted templates. The segmentation map is up- 
dated, assigning each pixel to the child with the largest 
contribution. 

Using the peaks previously found in the high resolution 
image, the deblending procedure is repeated for each u-K 
background subtracted, PSF matched image. Aperture 
photometry in 1.5" color apertures is then performed on 
the "deblended" children. Total fluxes are computed by 
dividing the parent total flux according to the aperture 
flux ratio of the children, thus conserving the total flux of 
the parent. The IRAC and MIPS photometry of the chil- 
dren is remeasured using the method described in ji4.3i 
but now using the new deblended segmentation map and 
the positions of the children. Due to the large PSF size 
of the GALEX data and the lack of a deblended w-band 
segmentation map, no attempt was made to deblend the 
NUV and FUV data. 

2186 and 3017 objects in the original catalogs were split 
into 5019 and 6803 children in AEGIS and COSMOS, 
respectively. New identification numbers of the children 
start at ID=30000. On average, the difference between 
the total isT-band flux of the brightest child relative to the 
parent increases with redshift, with a median difference 
of 0.1 mag at z ^ 0.5 and 0.3-0.4 mag at z > 1.5. 

Figure [S] shows an example of an object that has been 
deblended with the above algorithm. Although classi- 
fied as a single object by SExtractor, there is clearly 
a faint neighbor. The deblending technique recovers a 



well-described SED of the faint neighbor, with a notably 
different shape and photometric redshift compared to the 
brighter galaxy. It is worthwhile to note that, in case of 
a fainter neighboring object (e.g., child 2 in Figure [5]), 
the photometric redshift of the brighter object (child 1) 
is generally not significantly different from the parent, 
even though the children can have vastly different spec- 
tral shapes and photometric redshifts. 

4.6. Star Classification 

It is important to distinguish between galaxies and 
objects that are likely stars in the photometric catalog. 
Stellar SEDs have a characteristically steep rise and fall 
in their spectra, resulting in a tight stellar sequence in 
their u-J and J-K colors relative to the colors of galax- 
ies. This is shown in Figure [HI stars (red) form a distinc- 
tive sequence with blue J-K colors for a wide range of 
u-J colors as compared to galaxies (black). 

Stars are fiagged based on the u-Ji and Ji-K colors 
of all objects, where any objects satisfying the following 
limits (corresponding to the red line in Figure [9]) will be 
classified as a star: 

Ji - A' < -0.74 + 0.26(w- Ji) [u-Ji<3] 
Ji- K < -0.55 + 0.19(m - Ji) [3 < u - Ji < <]12) 
Ji -A' < -0.16 + 0. 10(u- Ji) [ii-Ji>4] 

The middle panel in Figure [3] shows the va lues for 
the st ars (red) and galaxies (black) when fit with iPicklesI 
(jl998[) stellar template library compared to the default 
set of EAZY galaxy templates. The vast majority of ob- 
jects that are classified as stars through their u-J and 
Ji-K colors are better fit with the stellar templates. The 
fraction of objects that are classified as stars relative to 
the total number of objects as a function of A"-band mag- 
nitude is shown in the right panel in Figure IH] for both 
fields. Almost all objects brighter than K ^ 16 mag are 
stars. 

We have compared our star classification to the COS- 
MOS Hubble Space Telescope (HST) Advanced Camera 
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Fig. 9. — Objects arc classified as stars (red) from their u-Ji and Ji-K colors in the left panel. These same objects also have SEDs 
that tend to be better-fit by the stellar templates than the galaxy templates (middle panel). All objects in both catalogs are shown, with 
no quality cuts applied. The fraction of stars relative to the total number of objects as a function of K-hand magnitude for both fields is 
shown in the right panel. 



for Surveys (ACS) catalo£3, matching objects within a 
radius of 1". Wc find that > 90% of the objects classified 
as stars using the higher resolution /-band data also fall 
within our u-Ji and Ji-K color selection limit in Equa- 
tion [131 The star classification method seems robust, 
although we note that there may be a small fraction of 
objects (< 10% of stars and < 1% of galaxies) that are 
misclassified, especially close to the detection limit. 

4.7. Catalog Format 

The final if-selected catalog consists of 24739 and 
27520 objects detected by SExtractor in the AEGIS and 
COSMOS mosaics, with a total of 27572 and 31306 ob- 
jects in the deblcnded catalogs. Photometry in the FUV 
through the MIPS 24 ^m band is included, with addi- 
tional optical broad- and medium-band photometry in 
the COSMOS field. The catalog fluxes can be converted 
into total magnitudes as follows: 

TO = -2.5 log F + 25 (13) 

where F is the total fiux density of the object with a 
magnitude m normalized to a zero point of 25 in the AB 
system, corresponding to a flux density of 3.631 x lO"'^" 
erg s~^ llz~^ cm~^ (or 0.3631 fJ,Jy). To convert from 
total magnitudes to the color aperture magnitude, the 
fltrxes and errors should be multiplied by the ratio of the 
-fC-band aperture fiux (Kaper) to the total if -band flux. 
A description of the columns included in the Version 5.0 
catalogs can be found in Tabled) 

A standard selection of galaxies can be obtained by 
selecting use=l for either the original SExtractor catalog 
or the deblcnded catalog, which is equivalent to wmin>0.3 
and star_flag=0. An additional cut on the S/N ratio 
of the galaxies may be optimal. 

4.8. Completeness 

We estimate the (point-source) completeness of our 
catalogs as a function of magnitude by attempting to 
recover simulated sources within the K-hand noise- 
equalized detection images. The simulated point sources 
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Fig. 10. — K-hand completeness curves; the completeness is de- 
fined as the fraction of simulated point sources that were recov- 
ered as a function of the total ii'-band magnitude of the simulated 
source. The point sources were inserted at random locations across 
the entire image where the exposure map was greater than half the 
maximum exposure, with no masking (solid lines) and detected 
sources masked (dotted lines). The inset plot shows the raw num- 
ber counts of galaxies within the two NMBS fields shown with 
Poisson error bars, with no corrections for incompleteness. 

were generated from the if-band empirical PSF de- 
scribed in § 14.11 before smoothing. The empirical PSF is 
scaled to the desired flux level and inserted at random lo- 
cations within the entire image. We then run SExtractor 
using the same settings described in § 14.21 determining 
the fraction of recovered sources as a function of input 
if-band total magnitude. Figure [10] shows the result- 
ing completeness curves from this analysis. The com- 
pleteness curves qualitatively agree with the raw number 
counts of galaxies within the fields shown in the inset 
plot. 

A fraction of the simulated point sources will fall on 
or close to other objects and these blended objects may 
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TABLE 2 

Summary of Photometric Catalog Contents 



Column No. 


Column Title 


Description 


1 


id 


Object Indentifier, beginning from 1 


2,3 


X, y 


X- and y-image coordinates 


4,5 


ra, dec 


Right ascension and declination (J2000.0; decimal degrees) 


6,7 


Kaper, eKaper 


Color aperture i^-band flux and error 


8-82 (8-67)'' 


X, eX, wX 


Total flux, error, and weight in each filter 


83-106 


X, oX 


Total flux and error of the COSMOS optical medium-bands 


107 (68) 


ftot24um_uJy 


Total MIPS 24 fim flux (^tjy) 


108, 109 (69, 70) 


f24um_ujy, e24um_ujy 


Aperture flux within 6" and error of MIPS 24 ^m (/jJy) 


110 (71) 


w24um 


Weight of MIPS 24 fj.m 


111 (72) 


wmin 


Weight of u—K band with minimum exposure 


112 (73) 


wminjrac 


Weight of IRAC band with minimum exposure 


113 (74) 


z_spec 


Spectroscopic redshift, if available (matched within 1" radius) 


114 (75) 


starjiag 


Star flag (l=Star, 0=Galaxy) 


115, 116 (76, 77) 


ap_col, ap_tot 


Size of the color aperture and total aperture in arcscconds 


117 (78) 


totcor 


Aperture correction applied to the iC-band AUTO flux to convert to total flux 


118, 119 (79, 80) 


K_ellip, K_theta_J2000 


ii'-band ellipticity and position angle from SExtractor 


120 (81) 


K.R50 


i^-band half-light radius from SExtractor (pixels) 


121 (82) 


K_class_star 


i^-band SExtractor parameter measuring stellarity of object 


122, 123 (83, 84) 


KJlags, UH2_fiags 


i^-band detection flag and maximum detection flag from U—H2 from SExtractor 


124 (85) 


Near_Star 


Signifies object is near a bright star and photometry may be affected 


125-130 (86-91) 


X_contam 


Ratio of the IRAC/GALEX flux in neighboring objects removed to the "cleaned" flux 


131 (92) 


contam_fiag 


IRAC/GALEX contamination ratio >50% in any band (1=>50% in >1 band, 0=OK) 


132, 133 (93, 94) 


nchild, id_parent 


Number of children from dcblended catalog (-1 if child) and ID of parent (-1 if parent) 


134 (95) 


use 


Standard selection of galaxies for either the deblended or the original catalogs (use = 1) 



^ AEGIS column numbers are in parentheses, as the COSMOS catalog contains additional bands 



not be properly handled by SExtractor. We have also 
repeated the simulations by inserting the point sources 
in locations that do not overlap with the segmentation 
map, thereby "masking" the original detected sources 
(dotted lines in Figure [TU]) . We include the 90% and 
50% completeness limits for both the masked and un- 
masked simulations in Table [3l We note that about 10% 
of the objects are blended within the original catalogs 
produced by SExtractor, as determined by the additional 
dcblcnding algorithm described in § 14.51 When using the 
deblended catalogs, these completeness limits may there- 
fore be slightly conservative. As the NMBS A'-band fields 
are not overly crowded, the completeness limits for both 
simulations are similar. 

4.9. Number Counts 

Figure [11] shows the surface density of galaxies as a 
function of magnitude in the NMBS fields compared to 
other galaxy number count analyses drawn from the lit- 
erature. Objects that are classified as stars (see § 14. 6p 
have been excluded from this calculation. The number 
counts are calculated in 0.4 mag bins using objects within 
the image area with >30% of the total A'-band exposure 
time (~0.21 deg'^ in both fields). 

The number counts of the NMBS fields have been cor- 



TABLE 3 
Point-Source Completeness Limits 
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Fig. 11. — Average galaxy number counts from the two NMBS 
flelds, along with a compilation of results from the literature. Pois- 
son error bars are shown for the NMBS flelds; the other fields have 
comparable uncertainties. The NMBS flelds have been corrected 
for incompleteness, but only points where this correction is < 20% 
are shown. 
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Entire Image 
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■ All magnitudes are given in the AB system. 



reeled for incompleteness using the simulations without 
masking, as calculated in § 14.81 We note that the com- 
pleteness values only provide a simplistic correction to 
the surface density of objects, whereas a more sophisti- 
cated analysis would account for the difference between 
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measured and intrinsic fluxes of tlie syntlietic sources. 
Due to tlie potential biases in the photometry of the 
faintest sources, we only extend the average number 
counts to K ^ 22.7, where the completeness correction 
begins to become significant (Figure [T0|). 

The number counts of the NMBS fields agree with the 
large compilation of data sets, within the scatter due to 
field to field variations. We note that the COSMOS field 
appears to have a relatively high de nsity of obj e cts at 
K < 18 mag. We refer the reader to IWake et alJ ()2011l ) 
for a detailed analysis of the relative densities and galaxy 
clustering in the NMBS fields. 

4.10. Comparison to Other Photometric Catalogs 

We have compared our K- and i^T.s'-ba nd flux mea- 
surements to the COSM OS survey catalog (jCapak et al.l 
[20nl lllbert eTHIIMffll FI. The COSMOS catalog is I- 
selected and quotes fluxes within a 3" aperture with an 
aper ture correction provided t o convert to AUTO fluxes 
(see iCapak et all 120071 I2008D . The Ks-hand data in 
the COSMOS catalog is derived from an earlier CFHT 
Ks image taken during the 2004-2005 observing season, 
whereas we use the publicly available Ks image with 
photometry from the 2005-2007 observing seasons for 
this analysis. As the COSMOS Ks image used to make 
the publicly available COSMOS catalog has 1.5" seeing 
(compared to ^ 0.8" seeing in the image version used in 
the NMBS catalogs), we smooth our Ks image with a 
Gaussian kernel to 1.5" and compare the fluxes between 
the catalogs. 

Figure [12] shows the difference in if^-band magnitude 
between the smoothed NMBS Ks image and the magni- 
tude quoted within a 3" aperture in the COSMOS cat- 
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Fig. 12. — The difference in X-band magnitude between the 
NMBS Ks imago smoothed with a Gaussian to 1.5" seeing and the 
COSMOS catalog magnitude within 3" apertures (black points). 
The running median (red) indicates that there are no significant 
systematic zero point offsets between the two catalogs. 

|http : //Irsa. Ipac . caltech. edu/data/COSMOS/tables/ 
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alog (black points) for all objects. The results indicate 
that there are no significant systematic zero point offsets 
between the NMBS and COSMOS catalogs. 

We have also compared the GALEX and IRAC flux 
measurements to publicly available catalogs, flnding ex- 
cellent agreement. We exclude comparisons to catalogs 
where it is unclear how aperture corrections are applied 
as the interpretation is difficult, but note that there were 
not any large systematic offsets. The COSMOS GALEX 
NUV and FUV catalog fluxes within 3" apertures agree 
within ^ 0.05 magnitudes with the public catalogs avail- 
able through the Multimission Archive^, The AEGIS 
NMBS IRAC photometry within 3.5" apertures has a 
median difference of < 0.03 magnitudes for all chan- 
nels when compared to the EGS cataloeP^. Finally, we 
have also compared the MIPS total fluxes in COSMOS to 
the sCOSMOS G03 catalog, finding excellent agreement 
with a median offset of < 1%. In general, the photom- 
etry of the NMBS catalogs agrees well with all publicly 
available catalogs within the fields. 

5. PHOTOMETRIC REDSHIFTS 

Photometric rcdshifts and res t-frame colors were dc - 
rived using the EAZY codS (iBrammer et al.l 12008V 
that fits linear combinations of seven templates to the 
broad- and medium-band SEDs. A K-hand apparent 
magnitude prior and systematic errors due to template 
mismatch are taken into account. The optimized tem- 
plate set used for the EAZY photometric redshift code 
contains seven templates, which is large enough to span 
a broad range of galaxy colors, while minimizing color 
and redshift degener acies. The d efault template set is 

(2001); 



described in detail in IBrammer et al. (20081): this set in- 
cludes five templates generat ed based on the PEGASE 
popu lation synthesis models ()Fioc fc Rocca-Volmerangj 
Il999t) and calibrated with synthetic photometry from 
semi-analytic models, as well as an additional young and 
dusty template added to compensate for the lack of ex- 
tremely dusty galaxies in semi-analytic models. We in- 
clude an additional template for an old, red galaxy from 
the lMarasto'^ (|2005[ ) models, with a lKroupal |2001f ) IMF 
and solar metallicity for a stellar population that has an 
age of 12.6 Gyrs. The choice of template set is impor- 
tant for both the photom etric redshifts and res t-frame 
colors (see Appendix C of IWhitaker et al.ll2010l for log 
M0 > 11). We adopt z_peak as the photometric red- 
shift (except where otherwise noted), which finds discrete 
peaks in the redshift probability function and returns the 
peak with the largest integrated probability. 

Uncertainties in the photometric zero points of the nu- 
merous photometric bands can lead to systematic shifts 
in the observed colors, and thus the redshift estimates, 
of the sample. We compute offsets of the photometric 
zero points by fitting the EAZY templates to the full 
optical-NIR SEDs of objects with spectroscopic redshift 
measurements, with the fit redshift fixed to the spec- 
troscopic value. The offsets are computed itcratively by 
minimizing the de-redshifted fit residuals, which allows 
the separation of zero point errors from systematic ef- 
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Fig. 13. — Comparison be tween the NMBS photometric redshifts of the d eblended catalog an d spectroscopic redshifts from the DEEP- 
2 sampl e lIDavi s et al.l [2003 ) and a smaller sample of high redshift LBGs l lSteidel et al.l [20031 ) (left panel), and the zCOSMOS survey 
l|Li llv ct al. 2007) (right panel). A few preliminary spectroscopic redshifts from Bezanson et al. [in preparation) and Labbc {in preparation) 
arc included. For reference, the red dashed lines are for Zphot = Zspcc it 0.05(1 + Zspcc). The correspondence between the photometric and 
spectroscopic redshifts is generally very good. 
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Fig. 14. — The scatter in the photometric and spectroscopic red- 
shifts of the deblended catalogs is relatively uniform as a func- 
tion of the ii'-band magnitude, with a larger scatter amongst the 
faint, high redshift LBGs (b lue). The spectr oscopic redshifts are 
from the 2COSMOS survey (Lillv ct al.l l20071 ) (top panel) and the 
DEEP-2 sample (Davis ct al. 2003) with a small sample of high 
redshift LBGs ifStcidcl ct al.. ,2003.) (bottom panel). We also in- 
clude preliminary spectroscopic redshifts from Bezanson et al. [in 
preparation) and Labbe [in preparation). The red dashed lines are 
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facts associated with the choice of templates, as a given 
rest-frame wavelength is sampled by different observed 
bands at different redshifts. The zero point offsets are 
listed in Table S] and included in the catalog. 

While the majority of the bands require only small off- 
sets (< 0.02 mag), some of the optical band s require off- 
sets of as much as 0.3 mag (see also Table 1 in lllbert et al.l 



TABLE 4 
Zero Point Offsets 





AEGIS'' 


COSMOS 


u 


-0.22 


-0.23 


B 




-0.07 


G 


-0.01 


-0.02 


V 




0.12 


R 


-0.01 


-0.03 


Rp 




0.03 


I 


0.03 


0.00 


Ip 




-0.34 


Z 


0.05 


-0.01 


Zp 




0.11 


J 


-0.03 


-0.08 


Ji 


-0.00 


-0.01 


J2 


-0.00 


-0.01 


J3 


-0.01 


-0.02 


H 


-0.03 


-0.02 


Hi 


0.01 


-0.01 


H2 


-0.01 


-0.03 


K 


0.00 


0.00 


Ks 


0.04 


-0.05 


IAA27 




-0.11 


/A464 




-0.16 


7/1484 




-0.11 


/yl505 




-0.02 


/A527 




-0.06 


/A574 




-0.12 


/A624 




0.04 


IA679 




0.24 


IA709 




0.15 


IA738 




-0.04 


IA767 




0.07 


IA827 




-0.12 


''All 


zero point 


offsets are 



in magnitudes, 
ZPeazy = ZPno 



defined as 

linal + offset 



l2009f ). We have verified that this technique reliably re- 
covers artificial zero point errors of 0.02-0.3 mag added 
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photometric redshift 

Fig. 15. — The linear (top) and logarithmic (bottom) photometric 
redshift distributions of AEGIS (black) and C OSMOS (grey) using 
the EAZY z_mc redshift (sec definition from I WittmanlMoa ) . We 
use the standard selection of galaxies (use=l) with Ktot greater 
than the 50% completeness limits from the unmasked simulations 
(23.6 and 23.4 mag in AEGIS and COSMOS). The spikes coincide 
with known overdensities. There are 12,632 galaxies at ^ > 1.5 
shown here. 



to the observed photometry. The CFHTLS li-band zero 
point is uncertain and known to have issues; a direct 
comparison between the CF HT u-band and S PSS sug- 
gests a difference of 0.1 mag (|Erben et al.|[2009t) . The de- 
rived M-band zero point offset may be partly influenced 
by the choice of template set. As a result, it is diffi- 
cult to separate true zero point problems from template 
mismatch. Furthermore, the it-band zero point is not 
well constrained as it is not bracketed by other filters — 
except the GALEX filters, which have their own known 
zero p oint uncertainies up to ~ 0.3 mag (e.g.. lllbert et alJ 
[200l . 

We find the photometric redshifts in COSMOS to 
be in excellent agreement with the magnitude-limited 
/-band Tab < 22.5 spectrosc opic sample of re dshifts 
from the zCOSMOS survey (jLillv et al.l [2007I ). with 
(T^/(l-h2;)=0.008 for 1509 galaxies (right panel. Fig- 
ure [131), where az/{l+z) is the normalized median ab- 
solute deviation, or cnmad = 1-48 x MAD. These val- 
ues are quoted only for Zspcc < 1, as there are very 
few spectroscopic redshifts currently publicly available 
in the COSMOS field beyond this redshift. We also find 
excellent agreement between the photometric and spec- 
troscopic redshifts for a larger sample of 2525 objects 
at Zspoc < 1.5 in AEGIS selected to a i ?-band limit 
of R ab < 24.1 from the DEEP 2 survey (jDavis et al.l 
[200l with a,/{l+z)=0m7 (left panel, Figure[13l). Both 
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Fig. 16. — The distribution of apparent X-band magnitude as 
a function of the photometric redshift (z.peak) for all galaxies in 
the NMBS fields (top) and the number of galaxies as a function of 
their photometric redshift (bottom). 



fields have very few catastrophic failures, with only ~5% 



> 5(7 outliers, defined such that Izphot 



:|/(l+^s: 



< 5(Tnmad- Spectroscopic redshifts are available for 
117 Lyman brea k galaxies (LBGs) a t z ^3 within the 
AEGIS field from lSteidel etall (pOOl . for which we find 
(T^/(l-|-z)=0.044, with 10% > 5a outliers. These rela- 
tively high phot ometric redshift un certainties for LBGs 
are well known ()Reddv et alJ 120081 ) and caused by the 
fact that LBGs tend to be faint in the rest-frame op- 
tical (observer's NIR) and their spectra have relatively 
weak Balmer/4000A breaks. We also note that there 
is excellent agreement between the NIR medium-band 
photo metric redsh i fts an d the Gemini/GNIRS redshifts 
from iKriek et al.l (|2006l ). with a biweight scatter in 
kphot - Z!,„c.r.\/(l + z .r,p.r) of Only 0.010, albeit this is only 
for four galaxies (see Ivan Dokkum et al1[2009l ). The pho- 
tometric redshift accuracies range from az/{i + z) ~ 
1% for galaxies with stronger Balmer/4000A breaks to 
az/{l+z) ~ 5% for galaxies with less defined breaks from 
the sp ectroscopic redshift sample. Recently, iKriek et alJ 
(|2011| ) find that the photometric redshift errors must be 

< 2% given the shape of the observed Ha emission line, 
as built from composite SEDs of the medium-band pho- 
tometry itself. 

In an effort to fill in the void of spectroscopic red- 
shifts for A'-selected galaxies above z ^ 1, follow-up 
spectroscopy for several massive galaxies at z ^ 1.5-3.5 
have been persued. Preliminary spectroscopic redshifts 
are in excellent agreement with the NMBS photometric 
redshifts (see Bezanson et al. in preparation and Labbe 
et al. in preparation) . We include the currently available 
preliminary spectroscopic redshifts in the right panel of 
Figure m 

In Figure [T31 we see that the scatter in the 
photometric-spectroscopic redshift comparison is rela- 
tively uniform as a function of the i^T-band magnitude. 
There is a larger scatter amongst the high redshift LBGs 
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Fig. 17. — The rest-frame UVJ diagram for NMBS galaxies with S/N > 8 in the X-band out to a redshift of 3. 5. The grey scale represents 
the density of points, with the lines indicating the separation between quiescent "red sequence" galaxies and star forming galaxies (both 
blue and red). The quiescent galaxies can be tr aced to the highest red shift interval 2.5 < z < 3.5. Apparently, the first quiescent galaxies 
stopped forming stars by that redshift (see also lMarchesini et al.|[2blOI ). 



(blue), which is not surprising given the spectral shapes 
and faint A"-band magnitudes of these galaxies. 

The photometric redshift distributions of both fields 
are shown in Figure 1151 using the EAZY zjnc redshift. 
This redshift value is drawn randomly from the redshift 
probability distribution, where the distribution of these 
redshifts for a given sample of objects very closely fol- 
lows the summed proba bility distribution of those same 
objects (jWittmanl 120091 ). The fields have independent 
spikes in redshift space, indicatin g the presence of known 
overdensitics (see, c.g.. l01sen et al.. ,2007.). The overden- 
sities in COSMOS at z < 1 are more prominent due 
to the addition of the twelve medium-band optical filters 
and may also be the partly due to the relative variance in 
the size of the structures from field to field. Additionally, 
we show the distribution of apparent K-hand magnitude 
as a function of photometric redshift (now z_peak, as 
throughout the rest of the paper) of all galaxies in the 
NMBS fields in Figure [H 

6. EVIDENCE FOR QUIESCENT GALAXIES TO 

Z - 3 

Galaxies within the local universe lie in two distinct 
classes: actively star-forming galaxies and "quiescent" 
galaxies with evolved stellar populations and little ongo- 
ing star-formation. These popu lations show strong bi- 
modal behavior in their colors (jKauffmann et all 120031 : 
iBaldrv et al.l 1200^ . where the red, quiescent galaxies 
have a more luminous characteristic magnitude and a 
shallower faint-end slope of the luminosity and mass 
functions compared to those of blue galaxies. 

Quiescent galaxies fall into a distinct "red sequence" 
in color- magnitude diagrams at low redshift. How- 
ever, at high redshifts dusty, star-forming galaxies can 
have similarly red U-V colors, "contaminating" the 



red sequence. Including the V-J rest-frame color 
can help separate out the dusty, star - forming galax- 
ics ("Labbc et al."2007^' wirvts et al.ll2007t iWilhams et al.l 
12009: Ilbcrt et al. 2009). At fixed U-V, passively evolv- 
ing galaxies have blue V-J colors, whereas dusty galaxies 
have red V-J colors. The location of the red clump of 
quiescent galaxies in the UVJ plane shifts towards bluer 
U-V colors towards higher redshift, while the V-J colors 
remain essentially unchanged. This change in rest-frame 
color with redshift reflects the intrinsic passive evolution 
of the stellar population. However, the uncertainties in 
the rest- frame colors will also increase with redshift, po- 
tentially washing out any intrinsic bimodality. These two 
effects have made detections of the bimodality difficult 
at high redshift. 

This bimodal color sequence has been shown to 
persist out to z ~ 2 both photometrically (e.g. , 
Labbe et al.|[2005HDaddi et al.|[2005HTavlor et al.l l2009al: 



Williams et al. r i2009t IWhitaker et all 12011" and spec 



trosc opically (e.g., iG iallong o et al.l 20051 iFranzetti et al.l 
120071 iCassata et alT 2008: K riek et al.ll2008a|) . A popula- 
tion of massive, quiescent galaxies at z ~ 2 is surprising. 
The ages of the stellar populations of these galaxies are a 
significant fraction of the age of the universe, which was 
only a few Gyr at these redshifts. The existence of mas- 
sive galaxies with strongly suppressed star formation at 
such high redshifts provides useful constraints on mod- 
els of galaxy formation (e.g., Baugh 2006). It is unclear 
at what point in cosmic time that this bimodal color se- 
quence emerges, due to limitations governed by the accu- 
racies of the photometric redshifts and rest-frame colors. 

The accurate redshifts and rest-frame colors of the 
NMBS are well suited for studies of the evolution of the 
bimodal color distribution of galaxies. In lBrammer et aP 
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Fig. 18. — Example quiescent SEDs selected based on their U-V and V-J rest-frame colors from the highest redshift bin in Figure (TT] 
sorted from bright to faint in /f-band magnitude. The best-fit EAZY template is shown in grey and the medium-band filters in red. All 
of these galaxies show a prominent Balmer/4000A break. 



(|2009f) . we selected 25,000 galaxies from the NMBS show- 
ing that the dust-corrected rest-frame U-V color dis- 
tribution is bimodal out to z ^ 2.5. Additionally, 
iWhitaker et al.l ()201ClD measured an increase in the in- 
trinsic scatter of the rest-frame U-V colors of massive, 
quiescent galaxies in the NMBS out to z ~ 2. This scat- 
ter in color arises from the spread in ages of the quiescent 
galaxies, where we see both quiescent red, old galaxies 
and quiescent blue, younger (post-starburst) galaxies to- 
ward higher redshift. Both of these studies were based on 
NMBS catalogs from the 2008A and 2008B data alone, 
whereas here we are using a deeper catalog containing the 
enti re survey data set . Also using the latest NMBS cata- 
log, iMarchesini et al.l ()2010ii showed that the most mas- 
sive galaxies (> 2.5x10^"'^ M0) at 3 < z < 4 are typically 
red and faint in the observer's optical with quenched star- 
formation, extending the evidence for quiescent galaxies 
beyond z ~ 3 with a sample of 14 very massive galaxies. 

Rest-frame U — V and V— J colors are derived by inte- 
grating the redshifted rest-frame filter bandpasses from 
the best-fit EAZ Y template. W e use the UBV band- 
passes defined bv lMafz ApellanizI (|2006( ) and the 2MASS 
J bandpass. This method for determining rest-frame 
fluxes is similar to that used by the COMBO-17 sur- 



vey described by (jWolf et all 120031 ). We find that in 
most cases this method pr ovides very similar results to 
the algorithm des c ribed b v Fvaii Dokkum fc PVa^ ([19961 ) 
and iTavlor et al.l ()2009af ). though the direct template 
fiuxes are more robust when narrower, closely-spaced fil- 
ters are used, as is the case here. A more detailed expla- 
na tion of the rest-frame c olor algorithm will be provided 
by iBrammer et al.l (|2011h . We note that the observed 
u-band is used in determining the rest-frame U-V colors 
for low redshift galaxies (z < 0.4). In order to limit bi- 
ases, we do not use the zero point offset which gives the 
best template fits in the rest-frame color determinations 
due to the large uncertainties. 

Figure fTTl shows the rest-frame U-V versus V-J dia- 
gram for the NMBS galaxies, selecting galaxies with a 
S/N> 8 in the i^-band. For the first time, we show the 
bimodal distribution of star-forming and quiescent galax- 
ies based on their rest-frame colors out to a remarkably 
high redshift of z ^ 3. A pronounced "quiescent clump" 
is visible at red U-V colors above the star-forming se- 
quence out to z = 2.5, with some evidence of a bimodal- 
ity at 2.5 < z < 3.5. This clump appears to shift to bluer 
colors at higher redshift. We note that the shape of the 
red clump of quiescent galaxies at < z < 0.5 extends 
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to bluer U-V colors, whereas the analogs of these galax- 
ies appear to absent at all other rcdshifts. This effect 
may be caused by low mass quiescent galaxies with lower 
metallicities that are below our A'-band detection limit at 
higher rcdshifts. We also note that the U~V, V~J plane 
would be populated differently in an optically-selected 
survey: there would be more faint, blue galaxies that fall 
in the lower left corner. The incompleteness of this flux- 
limited sample leads to an apparent increase in the num- 
ber of quiescent galaxies relative to blue, star-forming 
galaxies at the highest reds hifts. Using a mass-se lected 
sample of quiescent galaxies. IWhitaker et al.l (|2010( ) have 
shown that the fraction of quiescent galaxies falls below 
40% at z > 2. 

As stellar population parameters are depe ndent on the 
stellar population synthes i s models (e.g., seelWuvts et al.l 
[20071: IMuzzin et aLl 120091 : IConrov et al.l 120091 1201(iir ^ 
adopt the mo d el- independent selectio n criteria of quies- 
cent galaxies of iWilliams et al.1 (|2009f ) based on the rest- 
frame U-V a nd V-J colors. We hav e modified the selec- 
tion limits of IWihiams et al.l (|2009[) at z > 1, with the 
following adopted diagonal selection criteria, 

([/- F) > 0.88 X (F- J) +0.69 [z<0.5] 
{U-V)> 0.88 x{V-J)+ 0.59 [z > 0.5 ] (14) 

Additionally, we modify the limits in U-V and V-J 
such that, 

([/ - y) > 1.3, (y -J) < 1.6 [0.0 < z < 1.5] 
iU~V)> 1.3, (y -J) < 1.5 [1.5 < z < 2.0] (15) 
{U-V)> 1.2, (F - J) < 1.4 [2.0 < z < 3.5] 

where these constraints on the U-V and V-J colors pre- 
vent contamination from unobscured and dusty starform- 
ing galaxies, respectively. The clump of quiescent galax- 
ies is separated by only ^0.3 mag in U-V color from the 
blue, star-forming sequence at z ~ 3. We note that this 
offset is larger than the typical uncertainties in the rest- 
frame colors at z ~ 3 of < 0.1 mag. We speculate that if 
we were able to resolve enough galaxies at z ~ 3.5-4, we 
might begin to see the epoch where this color bimodality 
breaks down. 

In Figure [THl we demonstrate the quality of the SEDs 
of quiescent galaxies in the highest redshift bin of Fig- 
ure [ITl The K-hand S/N of these z ^ 3 galaxies ranges 
from S/N^ 35 for the brightest galaxy {K ~ 22 mag) 
down to S /N~ 9 for the faintest galaxies {K ~ 24 mag) . 
The majority of the galaxies shown in Figure [18] have 
little to no rest-frame UV emission, indicating little on- 
going star formation. We note that the rest-frame optical 
and NIR colors are dominated by an old stellar popula- 
tion, and that a comparatively low level of star formation 
is possible even on the "quiescent" side of the bimodal 
population. 

7. IMPROVEMENTS ENABLED BY THE NIR 
MEDIUM-BAND FILTERS 

Throughout this paper, we have emphasized the su- 
perior quality of the NMBS photometric redshifts and 
consequently, rest-frame colors. Here, we provide a more 
comprehensive look at the improvements provided by 
medium-band filters compared to the standard broad- 



band filters commonly used in high redshift photometric 
surveys. 

In § [6l we isolated a sample of quiescent galaxies out to 
z = 3.5 from the NMBS using the rest-frame U-V and 
V-J colors. Although these galaxies span a large range 
of S/N, their spectral shapes will unambiguously contain 
strong Balmer/4000A breaks. This sample of galaxies 
is therefore ideal to characterize the power of the NIR 
medium-band filters relative to the canonical broadband 
filters. 

We have generated broadband catalogs in both fields 
with the same depth as the medium-band images, as 
well broadband catalogs with the deeper WIRDs 
J and H photometry. The WIRDS J- and iJ-band im- 
ages have depths that are 1.2 and 2.3 mag deeper than 
the NMBS medium-band images in AEGIS (0.1 and 0.6 
mag deeper in COSMOS), enabling us to test how depth 
affects the accuracies of the photometric redshifts. The 
images with the same depth as the NMBS survey were 
generated by combining the J2 and J3 medium-band im- 
ages to create a J-band image (J = (J2 -|- Jz)/2), and 
the Hi and H2 medium-band images for an iJ-band im- 
age {H = (^1"+ i?2)/2) for both fields. The broadband 
catalogs contain FUV/NUV, ugrizJHK and IRAC pho- 
tometry in both fields and we compare these to medium- 
band catalogs with FUV/NUV, ugrizJiJ2J^HiH2K and 
IRAC photometry. 

The photometric redshifts are fit for the medium-band 
and broadband catalogs using the same parameter set- 
tings and templates as described in § [SJ Figure [TO] com- 
pares the median 68% confidence interval, (cr/(l -|-z)), of 
the photometric redshifts for all quiescent galaxies in the 
medium-band and broadband catalogs. The confidence 
intervals of photometric redshifts are appropriate given 
the true scatter in Figure [13] (at least for z < 1). It is 
therefore reasonable to believe them, and hence a com- 
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Fig. 19. — The median 68% confidence interval, {a/{l + 2:)), of 
the photometric redshifts of quiescent galaxies selected by the UVJ 
rest-frame color selection technique. The confidence intervals are 
calculated for medium-band catalogs and broadband catalogs that 
replace the Ji , J2 , J3 , Hi and H2 medium-band filters with J and 
H filters. 
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Fig. 20. — Examples of quiescent galaxies ranging from z = 0.27 
to z = 2.91, showing the rest-frame wavelength range the NIR 
medium-bandwidth filters span (red) as a function of rcdshift. The 
best-fit EAZY models are smoothed with a Gaussian with a width 
of 0.15/im/(l -I- z), the resolution of the medium-band filters. The 
medium-band filters begin to sample the Balmer/4000A region at 
z > 1.5. 



parison of confidence intervals for the larger sample is 
meaningful. We limit the sample to a completeness level 
of >90%, with Ktot < 22.2 mag. We note that the COS- 
MOS field includes roughly twice as many optical filters 
than the AEGIS field, including the twelve Subaru optical 
medium-band filters. The additional deep optical broad- 
and medium-bands improve the confidence intervals of 
the COSMOS NMBS photometric rcdshifts at z < 1 by 
a factor of two, with little effect at z > 1. For simplicity, 
we restrict this analysis to only those filters included in 
both fields where we can strictly test the differences be- 
tween a standard broadband catalog and the addition of 
the NIR medium-band filters. 

The median la confidence intervals of the broadband 
photometric redshifts of quiescent galaxies are up to a 
factor of two times greater than those calculated with 
the NIR medium-band filters. We note that the for- 
mal errors on the photometric redshifts will scale with 
the S/N of the detected objects, where the median la 
confidence intervals will be smaller for brighter samples 
within the NMBS catalogs. For example, the typical la 
confidence interval for a galaxy with K'tot=18-20 is 0.02, 
whereas this value increases to ^ 0.04 for A'tot=22. The 
distribution oi a/{l + z) for the quiescent sample within 
each rcdshift bin is roughly Gaussian; the inset panels in 



0.08 
0.07 
0.06 

+ 

^ 0.05 
b 

c 

•-6 0.04 

0) 

E 

0.03 

0.02 

0.06 
0.05 
? 0.04 
0.03 



Quiescent 


1 ' ' ' 1 


All Galaxies 






* Broodbond 






■ Deeper Broadband. 






• Mediumbond 




Sta^-fo^ming^ 
1 . . . 1 ^ 






I I I I I I I I I I I I I I I I I I I I I I I I 
Galaxies with Spectroscopic Redshifts 
^ Actuoi error 
• EAZY estimate of error 



I I I I I I I I I 




1.5 2.0 
Redshift 

Fig. 21. — The median 68% confidence interval of the photo- 
metric redshifts of quiescent (red) and star-forming (blue) galaxies 
selected by the UVJ rest-frame color selection technique. The 
confidence intervals are calculated for both the medium-band and 
broadband catalogs. The top right inset panels show the difference 
between the average best-fit EAZY templates of the quiescent (left) 
and star-forming (right) galaxies, smoothed with a Gaussian to the 
resolution of the medium-band filters. The bottom panel shows the 
median confidence interval of only those galaxies with spectroscopic 
redshifts compared to the actual errors (stars) as determined from 
the photometric-spectroscopic comparison in [js] with the LBGs 
indicated separately as a diamond. The best-fit EAZY template 
for the sample of LBGs at Zspcc = 2.5 is bluer than the average 
star- forming galaxy, with a weak Balmcr/4000A break (sec bottom 
right inset panel). The error for both quiescent and star- forming 
galaxies is up to a factor of two times smaller when using medium- 
band NIR data as compared to deep broadband data. 



Figure [1^1 show the best-fit Gaussian distributions. The 
width of the distributions range from 0.01-0.02 for the 
medium-band catalogs, compared to a factor of two in- 
crease for the broadband catalogs. 

The difference between the medium-band and broad- 
band catalogs becomes more pronounced for those qui- 
escent galaxies with z > 1.5, where the NIR medium- 
band filters start to sample the rest- frame Balmer/4000A 
break region. Figure [201 shows examples of quiescent 
galaxies ranging from z = 0.27 to z = 2.91 with the best- 
fit EAZY models smoothed with a Gaussian of width 
0.15//m/(l -|- z), the width of the medium-band filters. 
The rest-frame wavelength range that is sampled by the 
NIR medium-bandwidth filters (shown in red) shifts to 
shorter wavelengths towards higher redshift, beginning 
to sample the Balmer/4000A region at z > 1.5. 

Figure 1211 compares the median confidence intervals 
for the photometric redshifts of quiescent (red) and star- 
forming (blue) galaxies in the medium-band and broad- 
band catalogs, as selected from the UVJ diagram in 
Figure [ITl The difference between the broadband cat- 
alog and the addition of the NIR medium-band filters is 



Photometric Catalogs and Rcdshifts from the NMBS 



21 



most important for quiescent galaxies at z > 1.5. Deeper 
broadband photometry (squares in Figure [21]) improves 
the errors at z > 1.5, but only to a certain degree. The 
formal errors will decrease as the error bars decrease with 
deeper photometry but they seem to reach a floor around 
(t/(1 + z) ^ 0.05, at which point increased spectral reso- 
lution is vital. 

Quiescent and star-forming galaxies have different 
spectral shapes: the only prominent spectral fea- 
ture in the spectrum of the quiescent galaxies is the 
Balmer/4000A break whereas the star- forming galaxies 
typically have rest-frame UV emission with strong Ly- 
man breaks. The top inset panels in Figure \T[\ show the 
spectral shape of the average best-fit EAZY template for 
the quiescent and star-forming galaxies, smoothed with 
a Gaussian to the resolution of the medium-band filters. 
We note that the decrease in the width of the confidence 
interval at z = 3 for the blue, star-forming galaxies for 
the broadband catalog is the result of the Lyman break 
shifting into the optical filter passbands. 

The median confidence interval for the blue, star- 
forming galaxies are up to a factor of ^ 1.5 smaller than 
for the red, quiescent galaxies for the catalog that in- 
cludes the medium-band filters. The medium-band filters 
provide more similar photometric redshift uncertainties 
for the quiescent and star-forming samples (i.e. they re- 
duce the huge divergence seen at z > 2 for standard 
broadband filters), which can be important for cluster- 
ing analyses. We note that the average blue, star-forming 
galaxy still has a well-developed Balmer/4000A break, 
which is why the medium-band filters help even for those 
galaxies. Additionally, the deep optical broadband data 
will help constrain the photometric rcdshifts for those 
g alaxies with rest-frame UV emission. 

Ivan Dokkum et ahl (j2006D showed that the "median 
galaxy" at 2 < z < 3 with mass > 10^^ M0 is faint in 
the observer's optical, red in the observed NIR and has 
a rest- frame UV spectrum that is relatively fiat in F\. 
Therefore, the typical massive high redshift galaxy can 
either have a spectral shape similar to quiescent galax- 
ies or that of a red, dusty galaxy. For both cases. Fig- 
ure [21] demonstrates that the NIR medium-bandwidth 
filters are necessary to improve the confidence intervals 
of the photometric rcdshifts of the typical massive high 
redshift galaxy. In fact, the photometric rcdshifts arc 
more accurate when including the NIR medium-band fil- 
ters regardless of the spectral shape. 

In addition to comparing photometric redshift con- 
fidence intervals, we can also look at how the frac- 
tion of catastrophic outliers and the overall accuracy 
of those galaxies with spectroscopic rcdshifts changes in 
both fields when using the broadband catalogs versus 
the medium-band catalogs. We find that the fraction 
of catastrophic outliers increases by a factor of two in 
COSMOS when using only broadband filters with similar 
depth to the NMBS survey. By increasing the depth of 
the J- and iJ-band images by 1-2 mag, as is the case with 
the AEGIS WIRDS data, the fraction of catastrophic 
outliers is the same as that resulting from the medium- 
band photometry. The scatter, crz/(l + z), increases by 
15% and 30% in COSMOS and AEGIS for the deeper 
broadband catalogs, consistent with the EAZY estimate 
of the errors at these rcdshifts (sec Figure [ 



We note that the EAZY confidence intervals are not 
well calibrated. To test how representative the 68% 
confidence intervals of the photometric rcdshifts are of 
the actual errors, we consider the sub-sample of galax- 
ies with spectroscopic rcdshifts. In the bottom panel 
of Figure I2H we compare the actual errors from the 
photometric-spectroscopic redshift comparision in ^J5] to 
the median confidence intervals. At z < 1.5, the EAZY 
error tends to overestimate the actual error by ~ 20-40%. 
On the other hand, the sample of LBGs at z > 1.5 have 
actual errors that are almost a factor of three times larger 
than the EAZY errors (diamond in bottom panel of Fig- 
ure l^T]) . We note that the LBGs have spectral shapes dif- 
ferent from the majority of NMBS galaxies (see bottom 
right inset panel in Figure [2T]) . with weak Balmer/4000A 
breaks and very faint i^T-band magnitudes. The median 
i^-band magnitude of the sample of LBGs with spec- 
troscopic rcdshifts is 23.9 mag, compared to a median 
A'-band magnitude of 22.2 mag for the sample of galax- 
ies with Zspoc < 1.5 (see also Figure [14]). For brighter 
A'-band magnitude limits at z > 1.5, we currently have 
very limited spectroscopic information (see Sj5]). 

8. SUMMARY 

We present the relatively deep and wide medium-band 
NIR imaging of the NMBS survey, consisting of two 
- 0.21 deg2 fields within the AEGIS and COSMOS 
surveys. The observations were carried out as part of 
a NOAO Survey Program on the Mayall 4m telescope 
on Kitt Peak, using the NEWFIRM camera with five 
NIR medium bandwidth filters. The full details of the 
observations and data reduction are described in ^2.21 
and ^J2] respectively, including various optimizations to 
the images to improve the image quality. The astrom- 
etry is accurate within < 0.1" (0.3 pixels). The fi- 
nal combined images are constructed from individually 
registered, distortion-corrected, weighted averages of all 
frames; the flatness of the background and excellent qual- 
ity of the data is readily apparent in Figures [3] and [?] 
The reduced NMBS images and weight maps are pub- 
licly available through the NOAO archive and through 
the NMBS site. 

We combined the NMBS data with UV (GALEX), 
visible (CFHT and Subaru), NIR (CFHT) and mid-IR 
(Spitzer/IRAC) data to produce a public A'-selected pho- 
tometric catalog, which we make available to the com- 
munit-y]^. As about 10% of the objects detected by SEx- 
tractor are clearly blended objects (for example, see Fig- 
ure[5]), we use custom scripts to deblend the objects using 
a higher resolution A"-band image. The deblcnded cata- 
log is included in addition to the original list of objects 
identified using SExtractor. With the NMBS catalogs, 
we summarize our main findings as follows: 

1. The NMBS catalogs contain ~ 13,000 galaxies 
at z > 1.5 with accurate photometric rcdshifts 
and rest-frame colors. We demonstrate the 
excellent quality of the NM BS SEDs through 
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2. We find excellent agreement with available spectro- 
scopic redshifts, with crj(l + z) - 0.008 for 1509 
galaxies in COSMOS and cr^/(l + z) ^ 0.017 for 
2525 galaxies in AEGIS. It should be noted that 
the currently available spectroscopic samples are 
heavily weighted toward blue, low redshift galax- 
ies. Several programs for follow up spectroscopy of 
NMBS targets at z > 1.5 are currently underway. 

3. We show evidence of a clear bimodal color distri- 
bution between quiescent and star-forming galaxies 
persisting to z ^ 3, a higher redshift than has been 
probed before. 

4. The median 68% confidence intervals of the pho- 
tometric redshifts for both quiescent and star- 
forming galaxies is up to a factor of two times 
smaller when comparing medium-band catalogs to 
broadband catalogs due to the increased resolution 
of the five NIR medium-band filters. 

Follow up programs will continue to add to the legacy 
of the survey, as our understanding of these high redshift 
galaxies grows and evolves. Calibrations of the medium- 
band photometric redshifts at 1 < z < 3.5 will be possi- 



ble with the upcoming HST Treasury program, 3D-HST. 
This NIR spectroscopic survey will partially overlap with 
the NMBS fields. 
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